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TECHNICAL MEMORANDUM 1401 


DIRECTIONAL STABILITY OF TOWED AIRPLANES* 1 


By W. Sohne 


1. STATEMENT OF THE PROBLEM 


So fax, very careful investigations have been made regarding the 
flight properties, in particular the static and dynamic stability, of 
engine -propelled aircraft and of untowed gliders. Thereby the mechani- 
cal work of airplane control was considerably facilitated for the pilot j 
he obtained a certain assurance in case of operating errors, and he was 
enabled to do blind flying. Furthermore, it became possible for the 
airplane designer to predetermine with certainty the stability of an 
airplane he was about to develop, on the basis of theory, and to influ- 
ence that stability by constructive measures . 

In contrast, almost no investigations exist 2 regarding the stability 
of airplanes towed by a towline. With the freight gliders known in 
Germany, a glider train is usually unstable. A towed airplane does not 
automatically hold to the flight path behind the towing airplane; on the 
contrary, its pilot must continually and carefully observe the towing 
airplane and balance path deviations by control deflections. Evidence 
for the imperfect stability of a glider train is also provided by the 
fact that blind towed flight of long duration is impossible . If the 


*"Die Seitenstabilitat eines geschleppten Flugzeuges." Deutsches 
Ingenieur-Archiv, vol. 21, no. 4, 1953, PP* 24-5-265. Extract from the 
dissertation of the author (Braunschweig 1947)* Principal reviewer: 
Prof. Dr. H. Schlichting; coreviewer: Prof. Dr. H. Blenk. 

^As a pilot of freight gliders in the second world war, the author 
was able to collect detailed flight experiences which stimulated him to 
undertake the theoretical investigations here described. They were per- 
formed after the end of the war, with the permission of the British 
Ministry of Supply, in collaboration with the Institut fur Stromungs- 
mechanik der Technischen Hochschule Braunschweig. 

2 

The only report to which the author had access, by H. Solf (Die 
Seitenstabilitat eines geschleppten Flugzeuges, FB-Bericht 1306 , 1940), 
considered the directional stability of an airplane, the towing point 
of which was assumed to lie in a vertical straight line through the cen- 
ter of gravity, and the yawing oscillations of which were neglected. 

The results seem therefore not sufficiently applicable for the purposes 
of practice . 
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towing airplane is lost sight of, the pilot of the towed airplane is 
forced usually after a few seconds to release the towline. Even brief 
temporary inattentiveness may lead to badly disturbed flight positions. 

In flight with undisturbed view and in calm air, corrections for dis- 
turbances are performed partially intuitively and the pilot of the towed 
airplane becomes hardly conscious of them. In the case of squally 
weather, however, a high degree of attentiveness is necessary; for 
freight gliders, in addition, tiring manipulation of the controls is 
required. 

The disturbed motions occur in the following two forms : 

1. As a longitudinal oscillation (periodic deviations of the flight 
path upward and downward) . It can be observed if the pilots of gliders 
and freight gliders are inexperienced, and is, in part, caused by too 
violent elevator deflections. An experienced pilot has no difficulties 
in eliminating these disturbed motions, except in blind flight or in 
low-speed flight with far -rearward position of the center of gravity. 

2. As a self -excited lateral displacement oscillation which is 
coupled with a rolling and yawing. motion. This lateral displacement 
oscillation appears particularly pronounced in freight gliders in the 
case of high wing loading. It considerably impedes the work of the 
pilot and may cause the towline to break, or force its release in 
special flight attitudes, for instance, in low -speed or in climbing 
flight . 

Thus, the following report will aim at investigating the directional 
stability of the towed airplane and, particularly, at determining what 
parameters of the. flight attitude and what configuration properties 
affect the stability. The most important parameters of the flight atti- 
tude are the dynamic pressure, the aerodynamic coefficients of the 
flight attitude, and the climbing angle. Among the configuration prop- 
erties, the following exert the greatest influence on the stability: 
the tow-cable length, the tow-cable attachment point, the ratio of the 
wing loadings of the towing and the towed airplanes, the moments of 
inertia, and the wing dihedral of the towed airplane (which is a deci- 
sive factor in determining c^^j . In addition, the size and shape of 

the towed airplane vertical tail, the vertical-tail length, and the 
fuselage configuration are decisive factors in determining c n p and c yp> 

the yawing moment and side force due to sideslip, respectively. 

For an untowed airplane, the longitudinal and directional stability 
cannot be pushed so far that the airplane automatically holds to a cer- 
tain flight altitude and to a certain course; likewise, a towed airplane 
can hardly be required to hold automatically to the flight path behind 
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the towing airplane, especially for the changing conditions of take-off, 
climbing and cruising flight. Rather, one must strive to relieve the 
pilot of the towed airplane as much as possible of the work of control. 
For blind flying, one. must attempt to make orientation by the direction 
of the, towline possible - unless an instrument can be developed which 
would show the position of the towed airplane with respect to the towing 
airplane from cable direction and cable tension. 


2. INFLUENCE OF TOWED FLIGHT ON TOWING AIRPLANE 
AND FREIGHT GLIDER ON THE BASIS OF FLIGHT EXPERIENCES 


The effect of the towing connection on the towing airplane is entirely 
different from the effect on the towed airplane. The towing airplane is- 
more stable in towed than in untawed flight. It can fly more slowly, its 
control effectiveness is of longer duration, and disturbances are damped 
more rapidly. Consider, for comparison, a kite. The kite is stabilized 
by the drag and the mass of its long tail. The effect of a controlled 
towed airplane on the towing airplane is similar. 

In contrast to the towing airplane, the towed airplane becomes 
unstable in towed flight. Calculation and flight experiences, in accord- 
ance, show that this instability increases with increasing lift coeffi- 
cient, increasing wing loading, and increasing climbing angle. The ratio 
of the wing loadings, thus the ratio of the mutual lift coefficients, 
play therein sin important role. For reasons of efficiency, it would be 
desirable that both airplanes fly at lift coefficients corresponding to 
their optimum glide ratios . For obtaining equal behavior in the case of 
squally weather^, it would be desirable, in contrast, that the lift coef- 
ficients should be about equal. The experiences of practical flight 
operation, however, showed as the most favorable ratio of wing loadings 
of the towing airplane and of the towed airplane, with consideration of 
the flight properties in towed flight, approximately the value 1.6. For 
this ratio of wing loadings, the pilot of the towing airplane - in order 
to maintain soaring ability himself - must always fly so fast that the 
lift coefficient of the towed airplane does not exceed the value 0.8. 

The pilot of the towed airplane does not attain a flight attitude with 
higher c values where the control effectiveness is no longer sufficient 

for compensating the lateral displacement oscillation and the rolling and 
yawing motion connected with it. The desire for a low landing speed and 


5jf the wing loading of a towing airplane is three times as high as 
that of an empty towed airplane and the two c a values in cruising flight 
are 0.75 and- 0 - 25 , a vertical gust of 6.6 m/sec at 190 km/hr causes a modi- 
fication of the lift coefficient Ac Q of O. 5 J+ or 0.6 so that there results 

for the towing airplane a load factor in the gust of 1 . 72 ; for the towed 
airplane, in 'contrast, one of 3*2. 
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for a possibility of landing outside of airfields likewise leads to a 
lower wing loading of the towed airplane. 

The pilot of the towing airplane has only a slight influence on the 
flight-rpath correction of the towed airplane. He must only observe a few 
fixed rules; he must, for instance, stay above the permissible minimu m 
speed, the permissible minimum curve radius, and the permissible gl j dl ng 
angle at low speed. At higher c a values, the towed airplane has a more 

favorable gliding angle than the towing airplane and may in a given case, 
if the engines are throttled, even overtake the towing airplane. 


3 . SETTING UP OF THE EQUATION OF MOTION 
(a) Symbols 


Beside the coordinate systems defined in DIN L 100 , namely, the 
system fixed in aircraft (x^, y f , if), the system fixed relative to 

the flight path (x a , y a , z^) and the system fixed relative to the 

earth (x g , y g , z g ), in this report the experimental system (x, y, z) 

is used prominently. This system may be interpreted as a combination 
of the system fixed in aircraft and the system fixed relative to the 
flight path (fig. 1). 

Furthermore, G denotes the weight of the airplane, F the wing 
area, b = 2s the span, i x , i y , i z the radii of gyration, h the 

vertical distance between the towing point and the center of gravity 
(fig. 5), a the horizontal distance between the towing point and the 
center of gravity, l the length of the towline, and k s = l/s the 
relative cable length. 

As aerodynamic quantities, V signifies the undisturbed flight 
speed, v y the component of the disturbed speed in y-direction, o^, 

oiy., o> z the angular -velocity disturbances in the experimental system, 

a the angle of attack, p the angle of sideslip, 7 the flight-path 
angle (7 > 0 climbing flight). 

The following is generally valid: Motions are positive with the 
axes. Moments and forces are positive with the axes, except for c w 

and c & . These two are positive against the axes. 

In the following table, we compile the form coefficients and moment 
coefficients known from H. Scholz and F. WenkA 

^H. Scholz and F. Wenk, Jahrbuch l$Al der deutschen Luftfahrtfor- 
schung, p. I 159. 
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Lift 


Lateral force due 
to sideslip 


Lateral force due 
to yawing 


Lateral force due 
to rolling 


Yawing moment due 
to sideslip 


Lateral damping 


Yawing moment due 
to rolling 


Rolling moment due 
to sideslip 


Rolling moment due 
to yawing 


Damping in roll 


Cable force 

Cable rolling 
moment 


Cable yawing 
moment 


Aerodynamic form coefficients 
of directional stability 


Moment 

coefficients 


c 


a 


'IP 


dc q 

dp 


dc r 


-qz 


Ss “za/ V 


Ss£U X a/ V 


c np 




~nz 


Sscu fca/ V 



5c N 

Ss “xa/ V 



c zp - - 


Sc L 

dp 



Sc L 

dso^a/v 



5c L 

dst^a/V 



Aerodynamic form coefficients Cable moment 
of cable force and cable moments coefficients 


c s = Cw + 
c ls = c s~ 


c a tan a 

cos a + a sin nr. 
s 



„ _ _ -h sin a + a cos a 

ns s' i 
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(b) General Course of a Lateral Motion 

Assume a glider train flying steadily straight ahead. Due to a dis- 
turbance, the towed airplane gets into a position with respect to the 
towing airplane which would be characterized by the following angles 
(figs. 1 and 4): a angle between the undisturbed flight direction of the 

towing airplane and the towline, X angle between the projection of the 
x a -axis of the towed airplane onto the x g y g -plane against the x g -axis, 

P sideslip angle of the towed airplane against its flight path, p rela- 
tive wind-bank angle. The angle between the flight direction of the 
towing airplane and the x-axis of the towed airplane is = X tan 7+0. 

There originates an excited lateral displacement oscillation with the 
following relative motions with respect to the towing airplane: lateral 

motion v y of the center of gravity of the towed airplane in y direc- 
tion, rolling motion about the experimental longitudinal axis, 

yawing motion cog about the experimental vertical axis. These relative 

motions cause the following aerodynamic forces and moments, already known 
from the theory of directional stability of the untowed airplane: 

Lateral force due to sideslip Y^, lateral force due to yawing Y m7 , 
lateral force due to rolling Y,,^ 

Yawing moment due to sideslip Np, lateral damping N^, yawing 
moment due to rolling N my 


Rolling moment due to sideslip 
yawing L,,^ , damping in roll 



rolling moment due to 


In addition, there are the forces and moments stemming from the cable 
force : 

Lateral cable force Y s , cable yawing moment N s , cable rolling 
moment L g 


The lateral motion of the untowed airplane is known to be charac- 
terized by the four roots of the frequency equation, namely, by the real 
damping -in-roll root Ap which determines the quickly damped aperiodic 

rolling motion, a conjugate -complex pair of roots = Ug^ ± iw 2 ^ 

which determines the generally damped yawing oscillation (u^ damping 
factor, Wgj frequency factor), and the real spiral root A 4 which 

determines the generally very slow aperiodic spiral motion which may be 
negative, that is, stable, or positive, that is, unstable. 
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The lateral motion of the towed airplane, in contrast, is charac- 
terized by six roots. 5 it is characterized, as in the case of the untowed 
airplane, by the damping-in-roll root ^ and the pair of yawing- 

oscillation roots A 23 = U 23 - i w 23 * -^ n addition, one has here a further 

conjugate -complex pair of roots Aj^ = ± iw^ which characterizes 

the lateral center-of -gravity displacement oscillation and is almost 
always strongly excited (fig. 3 ), and a real root Ag which determines 

the static directional stability in towed flight. 

In what follows, we consider mainly the lateral displacement oscil- 
lation and the static stability in towed flight. 


(c) Assumptions for Setting Up the Equations of Motion 

For simplification of the calculation, the following assumptions 
are made: 

1. Since it is only to be investigated whether the disturbed motions 

are stable, indifferent or unstable, it suffices to assume the disturb- 
ances to be small and thus to put cos p « 1, sin (j. p., cos p *1, 

sin p « p, 

2. As in other stability considerations, one may investigate longi- 
tudinal and directional stability separately, on the basis of flight 
experiences for small disturbed motions, also in the case of the towed 
airplane . 

3 . The behavior of the airplane with fixed control surfaces is 
investigated . 

4. It is assumed that the towing airplane holds to the flight path 
undisturbed by the lateral motion of the towed airplane. For the pilot 
of the towing plane, even slight rotations about the vertical axis 
caused by the cable tension are clearly noticeable j he counteracts them 
almost automatically by hand, or they are corrected by the automatic 
pilot. Since, furthermore, mass and moments of inertia of the towing 
airplane are considerably larger than those of the towed airplane, this 
assumption is realized with good approximation. 

5- Finally, it is assumed that the moment coefficients increase 
linearily with the rotational speeds and and the lateral 

motion Vy$ this assumption is quite customary in directional -stability 

theory and has proved to be permissible. 

^This result of the following investigations is anticipated for the 
sake of better understanding. 
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(d) Setting Up of the Differential Equations 
of the Lateral Motion 

The general equations of motion for the towed airplane state t ha t 
the vectorial sum of aerodynamic and cable forces equals the forces of 
inertia and that the sum of aerodynamic and cable moments equals the 
angular -acceleration moments. For the investigation of the lateral 
motion, it suffices to consider the balances of moments about the x-axis 
and the z-axis, and the equilibrium of forces in the y-axis. In the 
following formulation of the differential equations, the acceleration 
terms and the terms of the aerodynamic forces and moments in the form 
given by H. Scholz and F. Wehk 6 are used. 

1. The lateral-force equation.- The mass acceleration is 
(G/g)Vjc cos y. The aerodynamic lateral force equals 
Ap + Yp + Y^ + . After introduction of the aerodynamic form 

coefficients this expression reads 


R ya air ~ ^ 


c a^ + P c q3 + — f“ ' c qx + 


sa bca, 
V 


scu za„ 


-qz 


or, because of a = o^, (u Jca = - ayp 


R 


ya air 


qF 


c a^ + 3 




sai x s<u z 



( 1 ) 


In addition, one has a cable lateral force. Due to its lateral motion, 
the towed airplane travels longer distances than the towing airplane which 
is assumed to fly steadily straight ahead. Thereby, velocity differences 
and mass accelerations originate which find expression in variable cable 
forces. In order to determine variation and order of magnitude of the 
velocity differences and accelerations, these quantities were investigated 
for an oscillation obtained by flying. The result of this investigation 
permits considering the speed of the towed airplane as constant. Taking 
a climbing or descending angle y into consideration, one then finds the 
cable force to be 


c W + A tan 7 ( 3F( c w + c a '^ an 7) 

S = - cos p = cos p 

cos ( 0 - X cos y) cos(cr - X cos y) 

After introduction of the cable -force coefficient c g = + c a tan y 

and because of cos p = 1, S is 

^See footno'te b on p. A. 
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S = ^ Fc s 

cos (a - X cos 7) 

and the sideward cable -force component is 


( 2 ) 


S ya = S sin (a - X cos 7) = qFc g tan (a - X cos 7) ( 3 ) 


Thus the lateral -force equation reads 


— VX cos 7 = qF 
g 


c„P + pc 


fa ) 

\ qP y « i x ) 


s<u x 


SO^ 


'qx 


"^qz + 


c s tan(q - x cos 7)J 


w 


In this equation, the variables X, p, p, ci^., o^, and a occur. 

The influence of cry is very slight. The quantities cr^, and 

the velocity v in y -direction are chosen as variables of the 

•7 0 O 

differential equation. Thus the variables X, p, p, and a are to 
be eliminated. For this purpose, there exist the following 
relationships : 


X cos 7 = 0^ - p y g = v yg = Vx cos 7 = vCo^ - p) 


■ V yp 

X cos 7 = —IS. 

V 


P -Ofc --2S 
v 


p - X sin 7 + X cos 7 sin p sin p = cr^ 

• • 

p = cjj + X cos 7(tan 7 - sin p sin p) = + X cos 7 tan 7 


P = ®>C 


yg 


tan 7 


Furthermore, there is (fig. 4 ) 
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-y g = i 

sin a + 

a sin (x cos 7 + p) 



= arc sin 

’_Zfi 

- — sin (x cos 7 + p) 





z 

Z 11 








\ 

= tan < 

arc 

sin 

- — - — sin (X cos 7 

+ 3 ) 

- X cos 7 > 




Z Z 








* 


tan (cr - X cos 7 ) = tan 


For this, one may put with good approximation 
tan (ff - X cos 7 ) = - 


y g / 

— + X cos 7 {] 

L + — ^ + P— 

_z l 

l) l 


(5) 


and thus 


— [tan (a - X cos 7 ) | = -(-££■ + + an— ^ 

dt L J \ z v ny 


(6) 


One eliminates the quantities X, p, p, and a by differentiating 
first the lateral-force equation (4). Thus one obtains 


* C af* + tan r ) + (“* - -r 110 


saw \ saw 

qp - i^q* + + 


S03 SO) 

3d, , 

y l x v 


'ZB + Ia + „ 1 

[l V z Z 


(7) 


Therein the term with dy is negligible . 

In equation ( 7 ), the directional-stability parameters, flight- 
mechanical time unit Tp = 2G/gpVF and relative airplane mass 

density p g = 2G/gpsF, are introduced, and s/V and z/V are eliminated 
by means of the following relationships: 


s = % 
V 


k = - 
s s 


l 

V 


T F k s 

^s 


l 

s 
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Then the equation ( 7 ) Is rearranged according to the variables m , a> , 

X Zr 

and v , and reads finally 

J O 


Tt 


Tx 


'-<^ c qz^ - “z( c q P - 'V J 


'<ix, 


Us 


c s|j - <% c qX^J - %Ca + 


V y§r } r , M/! c - c tan -k + c m. c M , Vi 
v F + V (°qP c a tan 7 + c s - «V c q X ~ J + ~ 


Z£o J±s_ . o (8) 
Vs 


2. The moment equations .- Rotary motions about the coordi na te 
axes are determined, for the untowed airplane, by Euler's equations for 
rigid bodies. For the towed airplane there occurs, in addition, a 
coupling of the rotary motions with the lateral motion of the airplane 
center of gravity due to the cable lateral force and its moments. Since 
Euler's equations are here referred not to the main inertia system but 
to the experimental system, they do not appear in the simple form known 
but in a general form, with the left sides? 


U x - J x^x + ( J z - J x) ta g 2aj(a z - a yj“ c 4ic(Jz “ J x ) + ^(Jy " J z ) 

4 

u z = J z<4z + (Jz - Jx) - «y|+%( J x - Jy) + ^(Jz - J x) tan g - 2 - a 

> 

(9) 

Therein J x , Jy, J z signify the moments of inertia about the experi- 
mental axes. Expressed by the main inertia moments, one obtains 
because of y = y^ 


J x = J ^ cos 2 a + J zf sin 2 a 


J z = sin 2 a + J Z £ cos 2 a 


J = J „ 

y yf 


(j v _ jJ tar^2a = = (j rf _ j^) silLJa 


2 


( 10 ) 


?H. Scholz, Jahrbuch 1940 der deutschen Luf tf ahrtf orschung , p. I 433. 
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The aerodynamic force moments L about the x-axis and N about the 
z-axis are composed of contributions stemming from the motions of side- 
slip, rolling, and yawing. They read, with use of the aerodynamic coef- 
ficients of these motions 


L = -qFs 


SOVo SOX,o 

P c Ip + ~T^lx + ~Y^lz 


N = -qFs 


saVn SdL„ 

3 c n P + “^nx + "Y^nz 


or, because of and cos p - my sin p 558 ac^ - oi^p 


L = -qFs 


N = -qFs 


pc 


ip " c lx" 


setyX sojjj. 


SUXj 


y + y T- z 


0 c «« - c 


n p nx v 


st V | + 


S(a z, 


nx 


-nz 


( 11 ) 


With the cable force ( 2 ) the lateral component of the cable force, 
referred to the experimental y-axis, becomes (fig. 4 ) 


Sy = S sin (a - X cos 7 - p) 


qFCs- 


sin (a - X cos 7 - p) 
cos (0 - X cos 7) 


For the cable rolling moment and cable yawing moment (fig. 5 ), there 
results therefore 


•^cable " +c L Fc s* 1 e 


s in (a - X cos 7 - p) 
cos (0 - X cos 7) 


N 

cable 


•r, sin (0 - X cos 7 - p) 
+ qFc a — 

s e cos (0 - X cos 7) 


Instead of using the quantities h e and a e , referred to the experi- 
mental system, the quantities may be referred to the system fixed in 
aircraft h and a with the transformation 
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h- e - h cos a + a sin a ag = - h sin a + a cos a 

and also the aerodynamic form coefficients of the cable moments c ls 
nnl c ns are introduced. Then the cable -moment equations read 


L cable = +l Fsc Zs 
N cable = +lFsc ns 


sin (a - x cos y - p) 
cos (cr - x cos 7) 

sin (cr - X cos 7 - p) 


( 12 ) 


cos (cr - x cos 7) 

With ( 9 ), (ll)* and (12) , Euler's equations now read 

+ ay(«T - J z ) 


J x«y + ( J z ~ J x) tan o 2 % - ay 


-“hc&z - J x )- 


-qFs 


P c Zp " c Zx' 


say\ say 


say 


'Zx 


J Zz ” c Zs 


sin (g - X cos 7 - p) 


cos (a - x cos 7) 


J z^z + ( J z - J x ) tan 2a oy - ay 


ay(J x - J y ) + ay(J z - J x ) 


tan 2a 


-qFs 


pc 


np " c nx" 


say\ say 


say 


sin (a - x cos 7 - p) 


v nx T mz ” ^ns' , . 

v V cos (a - X cos 7) 


(13) 


In these equations the directional-stability parameters T F and p 
are introduced as follows: 


Jx 

201/ 

qFs 

gp^Fs 

Jz 

= T//iz' 

qFs 

u \ s 
's \ / 




( 14 ) 
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Furthermore, the equations (13) are multiplied hy n s = and 

(l \2 A \2 

divided hy and ( ) , respectively. They then assume the form 


^t f 2 + 


J y " J z m p 

OJ — -T 2 

z j -t 1 
°x 


- 1 


= - P 


tan 2cti j 2,:, 


s 


F o^z “ i " “x 


/_ s 


'i \2 

r - 1 

^ x/ 


tan 2oti j 2 + 


Gzp \V ' c * x \v ^ 


2 F 


/_ s 


+ 


m rn c f ±-\ 2 - u C f S f Sln (P - X COS y - g) 

z F lz ^i x j s Zs y i x j COS (a - x COS 7) 


“z T F 2 + 


1 ‘ ~ 

Vi.- 


to-SOfr 2^ 




1 - 


'i Y \2 


m 2tan2a 
F 2 


= -<P 

' V 

2 


J x - J y T 2 


F “x - j T F + 


y(zh ■ 


a bc T F c nx^ z j + a Vz^F c nz|^ z j “ ^s c ns^ 


2 

s \ sin (a - X cos 7-3) 


cos (a - X cos 7) 


> ( 15 ) 


One simplifies these equations hy introducing the moment coefficients 
and the influence factors of form and position of the ellipse of inertia 


K l = 


tan 2a 


*n = 


/i x \ 2 

1 - -5 

UJ 


tan 2a K _ 1_ 


1 " “Z*n 

(16) 


then one differentiates them in order to he ahle to eliminate the 
expressions p, X, and a. Because of the value a (5) the expres- 
sion at the right in (15) becomes 
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sin (g - x cos y - g) 
cos (a - x cos 7) 


sin < arc sin s - sin (X cos 7 + 0 )- - X cos 7 - 0 


~ yg _ 

Z 

in r.2k 

z 


I J 

cos <arc sin --&■ - sin (X cos 7 + 6 ) 5 - -X cos y 
l l 


For this equation one may write with close approximation 


sin (q - X cos 7 - 0) _ 
cos (a - X cos 7) 


+ (x cos 7 + 0) (l + - 

i V z, 


With = X cos 7 + p there results finally for the differentiated 
quotient 


d j sin (a - X cos 7-0) 
dt cos (a - X cos 7) 


yg A a 

— s - ov. 1 + - 
Z k Z 


After elimination of the quantities 0, X, and a and introduction of 
the expression v/z = M- s T F k s > the differentiated equations (15) are 

rearranged and read finally 


“x T r + + + Vl T F 2 + “z T Frz ^VFj + I pl 1 s * I x°V T F + 1 b * 1 b ( 1 + 


— Vs - WFj + + K T f VA - ' **: 


“s t F 2 + <% T F k - a y t n T F^ + % + I Vty T F + n st 1 s ^ 1 + ^ + Vn T F 2 + ~ 


’ [^8 n P - Wf] + * K Tf “^ T F • < ^ L 7 J ^ r F - «»n* h 0 
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The braces in (17) are, because of the factor ay, almost equal zero and 
are neglected later on. 


In addition to the mass and angular acceleration terms ^~Vyg, 
and the aerodynamic forces and moments [~ v yg> j of the differential 

equations of the untowed airplane, one has for the towed airplane the 


cable forces and moments 


g 


)• 


3. Simplification of the differential equations and dimensionless 
representation of their variables .- The differential equations of the 
lateral motion contain a number of smaller quantities with ay and ay 

which, generally, may be neglected. Furthermore it is useful to repre- 
sent the variables o^., aig,, and Vyg in the following - dimensionless - 

manner 


“Vf * p “A * r 


Vf 2 ■ P 




‘Vv ? - i 


»>z T F 3 • * 


( T ys/ V ) l *s ' v 

’yg/ v ) |I s T F ' * 
(v /vV T 2 = v 

\ yg / ;s f 


Neglecting the correction factors and using a dimensionless representa- 
tion of the variables, one expands the three differential equations (8), 
( 17 ), and the lateral-force equation, by T F , and the last one, moreover, 

by n ; they then read 
s 


p + pl x + K^r + rl z + r|x e 


kp + pn^. + r + rn z + r^g 


+ ^s ' ^ + 


n p + njl + 


!) 1 - 


v2„ + vlJ-B- = 0 
P s k s 


vn D + vn — = 0 


^s _ 


> (18) 


-P c qx - P c a^s - ^ c qz " r P S ( c q|3 ~ c s^ + ^ + 
^( c qp " c a tan 7 + c s ) + vc s~- = 0 
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(e) The Frequency Equation 
With the statement 

p = Cp e Xt / T F r = C r e At / T F v = C v e At / T F 


one obtains from (l8) the following characteristic equation in the form 
of a determinant: 


M x 

+ \i B 

*0 + ^ + f) 

+ l sf B 

+ 

n p + “■(i + f) 

+ n a^ 

- c a |1 c 

~^ c qz ~ 

( C dP - C Sj) 

* 2 + *( C q0 ' C a tan r + 


= o (19) 


Computation of the determinant yields the following frequency equation 
of the sixth degree: 


AqTv 6 + ApA 5 + A 2 A 4 + AjA 5 + A J+ A 2 + A^h + Ag = 0 (20) 

This frequency equation is two degrees higher than that of the untowed 
airplane. The coefficients read,, after one has multiplied them by 
= l/(l - K-i* n ), in order to be able to put 



18 


NA.CA TM 1401 


A 1 “ K a [ ! x + "I”! + ”i + K n 1 z] + =qp ' =» 7 + =„ 

*2 * “civ's - Vz + ( c qp - c a *“» 7 + =s)(2 x - *1 ”x + n z - S’z) ' c qz( n p - K n z p) * 


C q*( 2 p - Vp) + 4 fl (“p ' Vp) + l*»(l + f)(“a ' Vs)} + °8j£ 


A 5 * "a j( c 
Vsf 
Vs 

- <= a ton 7 + c„) - V*) + 4a * Vp) - V( 1 *»P - Vp) ‘ c qx(»p“x * V*) -■ 

'*p - Vp) - c a W a ( n p - Vp) + 4g(l + f) f(Cqp - c a tan r + c B ) (n, - K n l B ) * - 

+ c a ( n p ' Vpj] + ^ [(“a - - (*s - Vs) c qx] + ^(tc - Vx + n z - V*)} 

a 4 - V.' 

( C qP - 

c a( n p l z ' l p°z) + + fj - c a ?](Vx - Z e n x ) + 

^ tan y + c,){l + &)(Vs - Vs) " “qxf 1 + f) ( l p n s - Vs) + ’ Vs) + 

^(Vz - Vx) + ^( l x n s - Vs) + j^j°e( n p - Vp) + ( c qp + c s)( n s - Vs) + c a(*s - «l"s)]j 

*5 * V.' 

+ ° B ) (Vs - Vs) + c s(Vp - Vp) + =qx( ! s n p - Vp) * c a(V B ' ^V)] + 


e .C->-v»)H)} 

a 6 - V" 3 ^ + l)( l s n p ' n » l p) 


The terms not underlined yield the coefficients of the frequency equa- 
tion of the untowed airplane. In addition, one has for the frequency 
equations of the towed airplane the underlined terms which stem from 
the cable forces and moments . The influence factors of the ellipse of 
inertia k ^ , K n are generally so small that they may be neglected. 

(See table I.) 


(f) Stability Conditions of the Frequency Equation 

For stability considerations, a numerical calculation of the indi- 
vidual roots is frequently not necessary; it suffices to determine the 
limits within which stability exists . The directional stability of a 
towed airplane is given, if all roots of the frequency equation have a 
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negative real part. For this, the condition^ is valid that all coeffi- 
cients and the following determinants (for the equation of the sixth 
degree) must be positive: 


D 1 = A 1 


% = 


*5 

0 


D 2 = 


Ao 

Ag 

a 4 

*6 


Ao 

Ag 


A 5 


0 

Ao 

Ag 

a 4 


D x = 


d 5 = 


A 1 Aq ° 

A 3 A-2 A-^ 

A]^ A^ 

Al Aq 0 0 

A 3 A l Aq 

Aij A 4 Aj Ag 

0 Ag A*j A 4 

0 0 0 Ag 


0 

0 

Al 

A 3 

A 5 


r ( 22 ) 


In the case of instability of the lateral displacement oscillation, only 


the determinant 
that stability may be attained, 
there results 


Dcj is negative. It must therefore become positive so 


If the determinant is computed. 


A A 6 - A.A, + . hh. - A 3 - k M Z > 0 

Aq A x A 4 - A 5 - A 3 A 2 + A 5 2 / /A 1 

Therein A]_A 4 - A^ is small (about 1. 5' percent) compared to 

-A^Ag + Aj 2 /a 4 . If the second-degree term in the numerator of this 
expression is expanded, the term A^ 2 Ag 2 is likewise small compared 

to the preceding terms ( 1.5 percent). If the small terms are neg- 
lected after expansion, there results as the condition for stability of 
the lateral displacement oscillation 

A 3 2 A£ - A 3 A 4 A 5 + AgAip 2 - 2A 1 A 5 a £ ^ 0 (23) 


8 H. Bilharz, 
19^1, P- 96. 
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An explicit representation and general discussion of this equation vith 
use of the expressions for the individual coefficients is not possible. 
However, if one wants to determine the stability limits for two variables, 
for instance, k s or c a against one obtains - since the individual 

coefficients contain linear terms with - a cubic equation from which 

one can calculate for every c & . and k g the pertaining l g. 

The stability condition for static stability reads 
A6 = K a (^s 5 / k s) c a( i s I1 p ~ n s z p) = 0 

Consequently, 


z s n p - Vp ^ 0 

If one replaces the moment coefficients by the aerodynamic form coeffi- 
cients and factors out various constant terms, this equation reads 

c np(jj + tan a ) * c Zp(“| tan a + lj >0 (24) 

If the towing point is located in the nose on a horizontal in front of 
the center of gravity ^h g =0, a e > 0^, there must be 

- c zp = +5c L / 5 P = 0 

This condition is to be fulfilled by a suitable dihedral. If the towing 
point is located above the center of gravity (a e =0, h e > o), there 
must be 


c np = " 5 c n/ 6 P ^ 0 

This condition is satisfied in the case of directional stability of the 
airplane . If the towing point lies at the center of gravity, static 
indifference results . 
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(g) Special Cases 

Since the work involved for the determination of the coefficients 
and for solution of the frequency equation of the sixth degree is very 
high, we investigated to see whether a few special cases with -a lower 
number of degrees of freedom could not serve as approximate solutions. 

In the calculation of the motion of the center of gravity of a 
towed airplane without consideration of the rolling and yawing motion 
(one degree of freedom), thus for = 0, there results a quadratic 

frequency equation, the conjugate-complex pair of roots of which charac- 
terizes the lateral displacement oscillation. 

In the calculation of the center of gravity and rolling motion of 
a towed airplane without consideration of the yawing motion (two degrees 
of freedom), thus for 0 ^= 0 , there results a frequency equation of the 

fourth degree, the real roots of which characterize the damping in roll 
and static stability in towed flight, and the conjugate -complex pair of 
roots of which characterizes the lateral displacement oscillation9. 

In the calculation of the center of gravity and yawing motion of a 
towed airplane without consideration of the rolling motion (two degrees 
of freedom), thus for = 0, there results a frequency equation of the 

fourth degree, the conjugate-complex pairs of roots of which characterize 
the yawing oscillation and the lateral displacement oscillation. 

In the calculation of the rolling and center-of -gravity motion of a 
towed airplane where the yawing motion is rigidly coupled with the center- 
of -gravity motion (two degrees of freedom), there results a frequency 
equation of the fourth degree, the real roots of which characterize the 
damping in roll and static stability and the conjugate-complex pair of 
roots of which characterizes the lateral displacement oscillation. 

A numerical comparison of the roots of the special cases with those 
of the general frequency equation showed that the lateral displacement 
oscillation and the static stability were not reproduced even approxi- 
mately correctly in a single case. The influence of the yawing motion 
as well as of the rolling motion on the entire directional stability is 
so great that none of those two may be neglected, without completely 
falsifying the results. Thus no approximate solutions may be found in 
this manner, and it remains necessary to seek the solution of the complete 
frequency equation of the sixth degree with all stated degrees of freedom. 

^This case is identical with the calculation performed by Solf (see 
footnote 2, p. 1) where coupling in a vertical straight line through the 
center of gravity was assumed and the yawing oscillation had been 
neglected. 
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if. DETERMINATION OF THE AERODYNAMIC AND 
CABLE FORCES AND MOMENTS 


The determination of the coefficients of the individual aerodynamic 
forces and moments is an important and difficult part of the calculation 
of directional stability. The coefficients dependent on the angle of 
sideslip, lateral force due to sideslip c^p, rolling moment due to side- 
slip c^p, and the directional stability c n p, are obtained by six- 

component measurements in the wind tunnel. The influence of modifications 
of the wing shape and of the dihedral, of shape and length of the fuselage 
and of the vertical-tail surfaces area, however, can be expressed numeri- 
cally in a simple form. 10 The coefficients dependent on the rolling and- 
yawing angular velocity, yawing moment due to rolling c^, damping in 

roll c Zx , lateral damping c nz , and rolling moment due to yawing c iz 

have, so far, been obtained almost exclusively by calculation.-^ The 
coefficients of the lateral force due to rolling c_„ and of the lat- 

eral force due to yawing c^ z are very small and are neglected below, 

as having insignificant effect on the results of the calculations. The 
cable-force and cable-moment coefficients may be obtained numerically 
in a simple form. The cable lateral force c s is a function of the 

aerodynamic efficiency and of the angle of climb of the towed airplane. 

The cable yawing moment c ns and the cable rolling moment c^ s are, 

moreover, dependent on the geometrical position of the towing point with 
respect to the center of gravity (fig. 5)* I n the case of the towing 
point lying in the nose the cable yawing moment is the decisive factor. 


H. Scharn, Systematisch-kritische Beispielrechnungen zur Seiten- 
stabilitat des Flugzeuges, Zentrale fur wissenschaftliches Bericht- 
swesen, UM-Bericht No. 2032, 1943; H. Schlichting, Aerodynamik der 
gegenseitigen Beeinflussung der Flugzeugteile , Monographie des M. 0. S. (A) 
Volkenrode, Reports and Translation No. 171 (176). Bericht 46/ 5 des 
Instituts fur Stromungsmechanik der Technischen Hoschschule Braunschweig. 

■^G. Mathias, Z. Flugt. Motor luf tsch. , 23> (1932), p. 224; H. Henn 
and H. Weiler, Untersuchung uber die Seitenbewegung von Flugzeugen, 
Vorabdruck zum Jahrbuch 194-3 der deutschen Luftfahrtforschung, 

Technische Berichte 11 (1944), No. 3* 



NACA TM 1401 


23 


5. SOLUTION OF THE FREQUENCY EQUATION 


In view of the complicated structure of the general frequency equa- 
tion it is difficult to make statements regarding the roots of the fre- 
quency equation, and thus on the lateral motion of the towed airplane , 
from the coefficients or from the aerodynamic and cable force parameters. 
One must therefore completely calculate numerous systematic examples in 
order to recognize the influence of a few variables of the flight atti- 
tude and of configuration parameters on the directional stability. First 
IP 

Graffe's method was used for the numerical solution of the frequency 
equation. This method permits calculation of the roots with arbitrary 
accuracy. However, it is laborious and does not offer a possibility of 
error control in the course of the calculation. For this reason approxi- 
mate solutions for determination of individual roots were looked for. 


(a) Approximate Solutions for the Roots 
of the Frequency Equation 

The numerical calculations were based on the aerodynamic and con- 
figuration data of the freight glider Go 24-2 (fig. 6 ). The coefficients 
and roots of the frequency equation are, for this airplane, for a flying 
weight of 5,500 kg, a lift coefficient c Q = 0.436, and a relative cable 

length k g = 6 , of the following magnitudes: A 1 = l8.7> A^ = 52.4, 

Aj = 316.1, A^ = 24.8, A^ = 74.7, Ag = 40.0, damping-in-roll 
root A^ = - 16 . 7 , yawing oscillation Agj = - 1.0 t 4.2i, lateral dis- 
placement oscillation Aj^ = +0.17 ± 0.56i, static stability Ag = -0.37* 

For other flying weights, lift coefficients, and cable lengths, too, as 
also for other model types, the orders of magnitude of the individual 
roots show proportions similar to those of the example mentioned. Thus 
they are, except for the small difference between A^ and Ag, numeri- 
cally very different. Hence it follows that the large roots A-j_ and 
Agj are determined mainly by the coefficients Aj_ to A^ whereas the 
small roots Aj^ and Ag, in contrast, are determined by the coeffi- 
cients Aj to Ag of the frequency equation. For the approximate 

solutions one makes use of this fact by breaking the frequency equation 
of the sixth degree ( 20 ) down into two cubic equations: 


• L 2 F. A. Willers, Methoden der praktischen Analysis. 
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A^ + A-jA 2 + AgA + A^ = 0 


A ^1+ 2 Ac Ag 
^ + —A + -2-A + — = 


0 


(25) 

(26) 


From the first cubic equation (25) one calculates the large damping-in- 
roll root Aj and the pair of roots of the yawing oscillation Ag^ > 

from the second cubic equation (26) the pair of roots of the lateral 
displacement oscillation A^ and the small static-stability root Ag . 

A first approximate solution for Ag is obtained from (26) with 

the aid of a monogram) it is improved according to Horner's scheme. 
(Since the roots A^ and Ag are of the same order of magnitude, the 

expression Ag = -Agy^A^ is too inaccurate as a first approximation. ) 

One can then calculate A]^ = u^ ± iw^ with the aid of the Vieta 
relations, applied for (26). One has 


2u 4 5 + ^g 


A 5 




(27) 


For Aj_, Scham-*-5 indicated the approximate solution 


A-z 

Ai = A l + — - - 2 - (28) 

A 1 A l 


This approximation is perfectly sufficient if Aj_ is considered by 
itself. However, since A-l will be required again for calculation of 
the other roots, it is advisable to correct A^ by means of the Horner 
scheme referred to (25)* 

Since all roots, except Agj = Ug^ ± iwgj, are known, one may employ, 
for determination of u 2 ^, the Vieta relation for the entire frequency 
equation (20) 

\ + 2u 25 + 2u 45 + \ = u 25 = — - Y " ~ - u 4 5 ( 2 9 ) 

^•^See footnote 10 on p. 22. 
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Furthermore , one has approximated according to the Vieta relations for 
equation ( 25 ) 

r 25 2 \ = " a 3 w 23 = / r 25 2 " u 23 2 

A comparison of the solutions calculated according to these approximate 
formulas and according to Graff e ' s method showed, for a number of examples 
for which the lift coefficient and the dihedral were varied, the following 
errors of the approximate solution: 



*1 

U 23 

w 23 

u 45 

v 45 

*6 

Maximum error, percent 

0.04 

-3-7 

1.2 

- 19.2 

2.9 

-4.2 

Mean error, percent 

'. 02 

-2.9 

1.1 

- 16.6 

1.7 

-1.4 


The agreement is satisfactorily exact except for the error in the 
excitation of the lateral displacement oscillation u^. For some 

investigations these approximate solutions are, therefore, not sufficient. 
One can considerably improve them by expanding for Ag Horner’s scheme 

by two stages of (20), and for A^ by three stages (to include the entire 

equation (20)). One may improve the values of the conjugate -complex pairs 
of roots r]^ 2 = -A^j/A^Ag and r 2 ^ 2 = -A^A-, by replacing Aj by a 

corrected A^ for which the following relation is valid: 



(A2A]^ - 0.97ApAcj 


0.895A 2 A 5 


+ 1 


■ 79-ApAej J 


(30) 


This equation was derived from the Graffe method. The corrected A^ 
is from 1 to 3-5 percent (on the average 1.5 percent) smaller than A^. 

The error still possible is now smaller than 0.2 percent. Finally, the 
real parts of the two conjugate -complex pairs of roots are to be calcu- 
lated with the aid of the Vieta relations for the entire frequency 
equation (20). 
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(b) Damping Times and Period of Oscillation 
of the Lateral Motions 


.The damping or excitation times and the period of oscillation of the 
lateral motions result from the roots of the frequency equation in the 
well-known manner. By damping time, we understand the time during which 
the motion is damped to the l/e-fold value, t^ = -Tyj is the damping 


time of the rolling motion, t 


23 


= -T-p/u^ the damping time, of the 


yawing oscillation, T^ = P er i° <i of oscillation of the 


yawing oscillation, t^ = the damping or excitation time of 

the lateral displacement oscillation (> 0 indicates excitation time), 
T^__ = 2nTpy/w^^ the period of oscillation of the lateral displacement 


oscillation. An important parameter for the judgment of an oscillation 
is, furthermore, T^ jt^ = 2itu/w, the number of oscil l ations to damp to 

the l/e-fold value. 


(c) Desirable Stability Properties in Towed Flight 

In order to facilitate the work of the pilot, it was desired to 
obtain the following stability properties: 1. The lateral displacement 

oscillation is to be damped or at least indifferent. 2. The period of 
oscillation is to be as large as possible. 3» Static directional sta- 
bility in towed flight is to be only small. As shown by the calculations, 
a high static directional stability in towed flight is always combined 
with a pronounced excitation of the lateral displacement oscillation. 
Furthermore, a high static stability decreases the control effectiveness. 
4. Finally, a maximum damping and maximum period of oscillation of the 
yawing oscillation are desired, just like for an untawed airplane. 


6. NUMERICAL CALCULATION OF THE DIRECTIONAL STABILITY 


The directional instability in towed flight, in the shape of an 
excited lateral displacement oscillation, appears in a particularly inten- 
sive form in the case of the airplane Go 242 1 ^ (fig. 6) with high wing 
loading, selected for the calculation of examples. Furthermore, detailed 

-*-H/ith this model, the author was able to collect detailed flight 
experiences for various climatic conditions and loading conditions. For 
unfavorable conditions, for instance, with high pay load and at high 
temperatures when the engine power fell off, or in the case of small 
airfields where the pilot was forced to lift off the ground too early, 
numerous false starts resulted which frequently led to failure of the 
airplane . 
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six-component measurements y exist for the Go 242, so that one may refer, 
not only to estimated or approximate calculated values, but also to meas- 
ured values of aerodynamic force coefficients . The following calculations 
1 £> 

were performed: Starting from a theoretical coupling at the center of 

gravity and an infinite length of cable, --we investigated the influence of 
the position of the towing point and of the cable length on the origin of 
the conjugate-complex pair of roots A]^ of the lateral displacement 

oscillation and of the real root Ag . Then the effects of the most 

important parameters of the flight attitude and of several configuration 
parameters on the directional stability were determined. Finally, the 
stability regions for variable flight mechanical and configuration param- 
eters were determined with the aid of Routh's stability criteria. 


(a) Influence of the Cable Length and of the Position of the Towing 

Point on the Origin of the Lateral Displacement Oscillation 

and on the Static Stability in Towed Flight 

The approximate formulas for the roots offer a certain clue for 
recognizing the influence of individual coefficients on certain roots. 

A survey of the structure of the coefficients themselves is obtained 
if one breaks them down into the determining factors of table I. 

An airplane towed with infinite cable length and the tawing point 
situated at the center of gravity, is distinguished from the untowed 
airplane only by a variation of the spiral root A^ which in towed 

flight, due to the cable force, becomes negative and therefore stable. 

The coefficients A^ and Ag and the roots of the lateral displacement 

oscillation are not yet in existence. If one lets the cable length 
become finite (fig. 8), there newly originates the coefficient A^, and 
from the root A^ a conjugate-complex pair of roots develops which 

characterizes the lateral displacement oscillation. The oscillation 
here is still damped. Therein 2u^^ = A^, that is, the damping of the 

oscillation is half as large as the real root for infinite cable length 
from which the oscillation arises. The period of oscillation is infi ni te 

^H. Gothert, Sechskomponentenmessungen an dem Windkanalmodell 
Go 242, Unpublished report No. 82228 of the Luf tfahrtforschungsantstalt 
Braunschweig, 1943; F. W. Scholkemeier , Weitere Sechskomponentenmessungen 
an dem Windkanalmodell Go 242, unpublished report No. 82239 of the Luft- 
f ahr tf or s chung sanstalt Braunschweig, 1943. 

1 ^The airplane data, parameters, coefficients, an d roots of the 
frequency equation compiled in the tables had to be omitted. Reference 
is made to the dissertation of the author, Braunschweig 1947. 
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for infinite cable length. In fact, it attains relatively large 
values for very large cable lengths, for instance, k s = 24, cable 
length Z = 290 m, period of oscillation is 123 sec. 

If one then lets the towing point deviate from the center of gravity 
horizontally up to the nose, one obtains the new coefficient A£ and the 
static -stability root Ag which increases at first rapidly and then more 
slowly (fig. 9)- The excitation increases about proportionally to it. 

The oscillation frequency grows in the same region to somewhat more than 
double its value. If one lets the towing point deviate from the center 
of gravity vertically up to the upper surface of the wing, Ag and Ag 

will likewise originate (fig. 9)- The excitation grows in proportion, 
and the oscillation frequency increases; however - corresponding to the 
smaller distances - not as much as in the case of the towing point 
travelling toward the nose. 


(b) Influence of the Cable Length, of the Flight Attitude, 
and of Configuration Parameters on the 
Stability in Towed Flight 

1. Variation of the towing-cable length 1 '? .- The relative cable 
length was varied from 3 to 12, corresponding to an absolute cable 
length Z of 36.7m to 147m. The relative cable lengths exceeding 

^Transition to cable length 0: In the transition to a very short 

towing cable, the cable angle a may become very large also for very 
small lateral-displacement and yawing oscillations. Then, however, the 
approximation for a (5) is no longer correct. If, finally, Z becomes 
zero, a free yawing oscillation about the center of gravity of the towed 
airplane is no longer possible. This case may be expressed by the two 
differential equations of the rolling motion and of the motion about the 
towing point on the towing airplane (coupled yawing and center-of -gravity 
path motion) . Similar to the case 1=0 is the Fieseler rigid-bar tow 
(H. Solf, Jahrbuch 1942 der deutschen Luftfahrtforschung, p. I 391 and 
K. Petrikat and E. Pleruschka, Jahrbuch 1942 der deutschen Luftfahrtfor- 
schung, p. I 404) and the three-point tow (W. Sohne, Die Seitenstabilitat 
eines geschleppten Flugzeuges. Der Dreiecksschlepp, Monography of the 
M.O.S. Volkenrode, 1946) where, due to the two cables acting, at a larger 
distance, on the wing nose, a free yawing oscillation is no longer possi- 
ble , either . 
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these amounts k s = 24, 48, 96 have only a theoretical significance and 

were investigated in order to observe the origin of the lateral displace- 
ment oscillation. With increasing cable length, the damping in roll Ap, 

and tlje yawing oscillation u^>j + iwgj remain almost unchanged; excita- 
tion and frequency of the lateral displacement oscillation + iwj^ 

and static stability in towed flight A£, in contrast, decrease (fig. 10). 
An increase of the length of the cable as a means for obtaining better 
stability, however, is advisable only up to a cable length of about 
k g = 6 ( 76 m). For reasons of practical flight operation, a short towline 

is desirable. A cable length of k s = 5 - 6 (60 - 7^ m ) seems therefore 

to be most favorable. 

2. Variation of the lift coefficient .- Here we calculated with six 
different c a values from O .303 to 1.0, in order to cover the entire 

lift range. The lift coefficients above 0.681 appear, in practical flight 
'operation, generally only in the unsteady accelerated flight after the 
take-off. With increasing lift coefficient, the excitation u^ and the 

frequency w^c of the lateral displacement oscillation as well as the 
static stability in towed flight A 5 increase . The period of oscilla- 
tion Tjf slightly increases, up to the lift coefficient 0 . 6 , and then 
*45 

again decreases slightly whereas the excitation time t.., decreases 

^45 

considerably percentagewise (fig. 11). The pilot of the towed airplane 
is forced to balance by control deflections the instability which 
increases with increasing lift-coefficient values while simultaneously 
the effectiveness of the control surfaces decreases in proportion 
with the decreasing dynamic pressure. This is a difficult and very 
tiring task, particularly in squally weather. It explains the high 
degree of dependence of the stability on the flight speed. The calcula- 
tion therefore confirms the flight experience according to which a 
directionally unstable airplane can no longer be controlled, above a 
certain lift coefficient (for Go 242, for instance, c a > 0.8), in towed 

flight, although it exhibits, for this lift coefficient, still very satis- 
factory properties in untowed flight. 

3. Variation of the flying weight .- Of the possible actual loading 
conditions, we selected as examples the empty flying weight with trimming 
ballast of 4,200 kg, the standard flying weight of 5,500 kg, and the full- 
loading flying weight of 6,800 kg. With increasing flying weight, the 
excitation and frequency of the lateral displacement oscillation and also 
the static stability increase very considerably. For the excitation time 
and period of oscillation as well as for the damping time of the static 
stability there results a slight decrease, due to the influence of the 
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flight -mechanical time unit Ty (fig. 12). With increasing weight, the 
c a values, however, increase; the control effectiveness, therefore 
decreases. -Furthermore, the reserve power of the towed airplane decreases, 
and the flight speed is therefore lower than in the case of smaller weight. 
The impediments resulting from the increase in lift coefficient appear 
therefore once more, in addition. The difficulties for the pilot in 
balancing the excitation increase, therefore, considerably more than is 
expressed in the amounts of the excitation times and period of oscillation. 

4. Variation of the angle of climb .- The angles of climb selected as 
examples 7 = 3.44° and 6 . 85 ° correspond, for G = 5,500 kg and 

c a = 0.681 to the climbing speeds 2.8 m/s and 5.6 m/s. With increasing 

climbing speed, the excitation time of the lateral displacement oscilla- 
tion and the static stability in towed flight decrease only slightly 
(fig. 13). The frequency of the lateral displacement oscillation, in 
contrast, increases considerably with increasing climbing angle, that is, 
the period of oscillation decreases which is a phenomenon that is likewise 
inconvenient for the pilot. Since climbing flight is, of course, always 
combined with a lower flight speed, here again the phenomena of the 
increasing excitation of the lateral displacement oscillation in low-speed 
flight and of the decreasing' period of oscillation in climbing flight are 
additive. The calculation, therefore, confirms in this case, too, the 
flight experiences which show precisely for climbing flight at low speed 
the maximum instability and therewith the greatest difficulties in con- 
trolling the airplane . 

5 . Variation of the rolling moment due to sideslip by variation in 
dihedral .- The dihedral was varied from + 6 . 7 ° to - 5 . 3° whereby the rolling 
moment due to sideslip changed from c^p = -0.26 to +0.09- With 

decreasing dihedral, the excitation u^ of the lateral displacement 

oscillation decreases more and more, and after the zero point has been 
passed, it turns into a rapidly increasing damping (figs. 14 and 15 ). 

The frequency of the oscillation shows a minimum for excitation zero and 
then increases again. The static stability decreases about proportionally 
with the excitation of the lateral displacement oscillation and makes the 
transition to the unstable region at approximately equal dihedral. From 
this remarkable dependence of the directional stability in towed flight 
on the rolling moment due to sideslip, one may draw the following conclu- 
sion: The optimum flight properties result in the case of coupling at 
the nose for a rolling moment due to sideslip c^ = 0 ; for this moment 

the excitation of the lateral displacement oscillation is very small, 
the period of oscillation shows a maximum, and the static directional 
stability in towed flight simultaneously has become zero. For a high- 
wing monoplane, the rolling moment due to sideslip disappears at a 
slightly negative dihedral (Go 242: v = -2.23°). It is true that the 
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airplane in untowed flight then shows a certain tendency toward spiral 
dive which is, however, unimportant since the untowed flight always lasts 
only a few minutes . The primary object of consideration for a towed air- 
plane always is towed flight. 

6 . Variation of the vertical-tail surface area . - The. vertical-tail 
surface area was increased from 2 X 3*85 to 2 x 5 . 39m^ in the manner 
indicated in figure 6 whereby the lateral force due to sideslip, the yawing 
moment due to sideslip, and the damping in yaw increased considerably. To 
increase the vertical-tail surface area is a customary means for improving 
the directional stability of an untowed airplane. In the case’ of the towed 
airplane, there results likewise a large increase in damping, but - it is 
true - also a slight increase in the frequency of the yawing oscillation. 
The excitation and frequency of the lateral displacement oscillation, in 
contras-t, decrease only slightly (fig. 16 ) ; therefore excitation time and 
period of oscillation show a slight increase. For this reason, an increase 
in vertical-tail surface area would therefore hardly be a profitable means 
for improvement of directional stability in towed flight. However, the 
pilot has to balance the self -exciting lateral displacement oscillation by 
suitable aileron and rudder operation. This will be easier the more 
effective (and thus larger) these tail surfaces are. Thus, an increase 
in vertical-tail surface area is, after all, wholly suitable for improving 
the directional stability of a towed airplane. 

7 • Variation of the damping in roll . - In the. case of a variation of 
the moment coefficient of the damping roll from l x = 15.12 to 18.48 (as 

may be caused, for instance, by a variation in the aspect ratio of the 
wing from A = 8 to I 3 .I), the damping-in roll root A^ increases pro- 
portionally to the moment coefficient of the damping-in-roll 2 X . The 

excitation and the frequency of the lateral displacement oscillation 
decrease slightly, the excitation time and period of oscillation there- 
fore slightly increase (fig. 17 ). According to the calculation, the 
towing properties of an airplane with high damping in roll, that is, with 
high aspect ratio, are therefore somewhat more favorable than those of an 
airplane with low damping in roll, that is, low aspect ratio. On the 
other side, there is, of course, the greater inertia with which an air- 
plane with higher aspect ratio reacts to aileron deflections. 

\ 

8 . Coupling at the lower side of the fuselage and variation in the 
position of the towing point in horizontal direction . - Indifference of 
the lateral displacement oscillation combined with static indifference 
may be obtained without variations in the configuration of the airplane 
when the towing point is placed on the lower side of the fuselage, ahead 
of the center of gravity. In the case of coupling vertically below the 
center of gravity at the lower side of the fuselage, there results static 
instability and a damped lateral displacement oscillation. If one then 
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lets the towing point travel forward to a position 0.90m ahead of the 
center of' gravity, the root of the static stability Ag passes into a 

stable region. Simultaneously, the damped lateral displacement oscilla- 
tion at this point becomes an excited oscillation (fig. 18) and the 
oscilla'tion frequency shows a minimum. By coupling below the center of 
gravity, a tail-heavy moment is created which must be balanced by ele- 
vator trimming. A coupling at the lower side of the fuselage 0.90m 
ahead of the center of gravity seems, therefore, to be considerably more 
favorable than the customary coupling at the nose. 


(c) Stability Limits for Various Cable Lengths, Weights, and Lift 
Coefficients Against the Rolling Moment Due to Sideslip 

The variation of c^p for G = 5>500 kg, c a = 0-4-36, k s = 6 

resulted, in the case of coupling at the nose, in a special dependence 
of static stability and of the excitation of the lateral displacement 
oscillation on c^o * It appeared therefore desirable to check this 

result for other weights, cable lengths, and lift coefficients as well. 
With the aid of Routh's stability criteria, we determined the stability 
limits for dynamic and static stability for these variables against the 
rolling moment due to sideslip. According to this, the limit of static 
stability lies at = 0. For all weights the limit of static stability 

varies, as the cable length is decreased from a very small negative value 
to a somewhat larger positive value of 7p . However, the difference is 

so slight (approximately l/4° dihedral) that for all practical purposes, 
one may state that in the case of variation of the dihedral, for cruising 
flight c a = 0.436 the limits of static and of dynamic stability coin- 
cide at Z p = 0 in such a manner that one can obtain either only static 

stability and dynamic instability or vice versa, or static and dynamic 
indifference. In contrast, a simultaneous variation of the lift coef- 
ficient and of the dihedral (fig. 19) shows that below a lift coefficient 
of c a = 0.42, a gradually widening region of static and dynamic stability 

results. Above c a = 0.42, there is a region of simultaneous static and 
dynamic instability. 


7. SUMMARY 


With the aid of the well-known directional-stability theory 1 ® (for 
untowed airplanes), the directional stability of an airplane towed by 


18 


See footnote 4 on p. 4 and footnote 11 on p. 22. 
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means of a cable was theoretically investigated according to the method 
of small oscillations, and the . frequency equation for it was derived. 
Whereas the frequency equation of the directional stability of an untowed 
airplane is an equation of the fourth degree, there results for the towed 
airplane an equation of the sixth- degree . The damping-in-roll root A-j_ 
and the pair of roots of the yawing oscillation A^ = u^^ ± iw^ of the 

untowed airplane remain valid almost without modification whereas the 
spiral root A 4 , in contrast, disappears. In its place there appears a 
new conjugate -complex pair of roots Aj^ = uj^ ± iw^ which charac- 
terizes the lateral displacement oscillation occurring in towed flight, 
and a real root Ag which represents the static directional stability 

in towed flight. The lateral displacement oscillation is generally 
excited. The limit of dynamic stability may be determined by the Hurwitz- 
Routh stability criteria. We investigated what discriminant is decisive 
for the stability limit of the lateral displacement oscillation, and 
found that it was possible to simplify the corresponding expression con- 
siderably. The investigation of special cases in which one or two degrees 
of freedom (rolling and yawing oscillation) had been neglected did not 
yield usable approximate solut 5 ons . 

Numerous examples were calculated for the airplane Go 242. For 
coupling at the center of gravity, if one passes from infinite cable 
length to a cable length which is finite, there originates a damped 
oscillation from the real root A^ . If the towing point deviates from 

the center of gravity in the horizontal direction forward or in the 
vertical direction upward, this oscillation becomes an excited oscilla- 
tion and its frequency increases . Simultaneously, there appears a real 
root which represents the static stability in towed flight. Then various 
parameters of flight attitude and configuration were varied. It was 
found that in the present model the lateral displacement oscillation is 
strongly excited. Its duration amounts to 19 to 32 seconds, according 
to the cable length, thus to about the "J- to 10-fold multiple of the 
yawing oscillations. The most noteworthy results of a variation of the 
individual parameters are: 

1. With increasing cable length, the frequency and excitation of the 
lateral displacement oscillation decrease. However, it is not worth- 
while to increase the cable length beyond k s = 6 (l = 76 m). 

2. With increasing lift coefficient, the excitation and frequency 
of the lateral displacement oscillation increase. 

3- The same occurs in the case of increasing pay load. 

4. With increasing angle of climb, only the frequency of the lateral 
displacement oscillation increases. 
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These results confirm the flight experiences and flight measurements 
and explain many difficulties in towed flight and accidents that have 
occurred. They are probably similar for other model types with a high 
wing loading. 

5. In the case of variation of the dihedral, the excitation of the 
lateral displacement oscillation reverts, for a vanishing rolling moment 
due to sideslip (v = - 2 . 23 °) to zero, the period of oscillation here 
attains a maximum. Simultaneously, the static stability in towed flight 
becomes equal to zero. With the aid of the Routh discriminant, the sta- 
bility limits against the rolling moment due to sideslip were determined 
also for other weights, cable lengths, and lift coefficients. There 
resulted below c & = 0.42 a region of simultaneous static and dynamic 

stability. For this structural design (namely = 0 to 0.2, 

v = - 2 ° to 2 . 3 °), the towing conditions are therefore the most favorable 
ones. An absolute stability in towed flight in the sense that the air- 
plane is controlled by the cable forces and moments in such a manner that 
the pilot ' s work is made completely superfluorous is unattainable . For 
Zp > 0 the lateral displacement oscillation is damped, but static insta- 
bility exists whereas for < - 0.2 static stability exists, the lat- 
eral displacement oscillation, however, is excited. 

6 . With increasing area of the vertical-tail surfaces, the frequency 
and excitation of the lateral displacement oscillation decrease only 
little j however, under certain circumstances such an increase may be 
advisable for improvement of the control effectiveness. 

7 . An increase in the moment coefficient of damping in roll l 

x 

yields a small decrease of excitation and frequency of the lateral dis- 
placement oscillation. 

8 . In the case of coupling at the lower side of the fuselage 0.9Qm 
ahead of the center of gravity, static and dynamic indifference may like- 
wise be attained without configuration changes in the airplane. This 
position of the towing point seems therefore to be more favorable t han 
the customary coupling at the nose. 


Translation by Mary L. Mahler 
National Advisory C ommi ttee 
for Aeronautics 


TABLE I— COEFFICIENTS OF THE STABILITY EQUATION , ARRANGED ACCORDING TO DETERMINING FACTORS WITH NUMERICAL 
EXAMPLE FOR G = 5,500 kg; k 8 = 6; c a = 0.436 
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A a - coefficient* for the untoved airplane 

AA, « additional tern, 1 teaming fran the cable force, for towing point at the center of gravity, but Independent of the cable length 

AAy « additional term, • teaming from the cable momenta, but independent of the cable length 1 ^ 

Mj * additional term, steaming from the cable force at the center of gravity due to finite cable length 

Mjy » additional tern, a teaming from cable moment* due to finite cable length 

19 In the case of coupling at the nose, for c a - 0.56, the towing point lies at the same height as the center of gravity. Hence, 
he and l B * 0. For other angles of attack, too, l B re ma i n s small compared to rig, and was neglected for the sake 0 / simplification 
of the calculation. 
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Figure 3.- Excited lateral displacement oscillation of a towed airplane 
at a flight speed of 206 km/h and a cable length of 80 m. Towing 
airplane He 111, freight glider Go 242. Period of oscillation 
~ 32 sec. 
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Figure 4.- Cable forces and angles on the towed airplane. KP 
coupling point on the towing airplane; SP center of gravity of 
the towed airplane; FP towing point on the towed airplane. 
(The NACA reviewer suggests that the labels x , R cable, 

a ^a 

Ry. cable should possibly be , S„ , and S respectively.) 
y a. y a y 
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Figure 5.- Coordinates of the towing point on the towed airplane. 

h > 0 towing point above the center of gravity; a > 0 towing point 
ahead of the center of gravity (in the figure above, h is negative), 
hg = h cos a + a sin a, slq = a cos a - h sin a . 




Figure 6.- General view of 




Figure 7.- Polar of the freight glider Go 242. 
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Figure 8.- Origination of the roots of the lateral displacement oscil- 
lation for change of the cable length from infinity and for coupling 
at the center of gravity. G = 5500 kg, c a = 0.436. 
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Figure 9.- Roots of the lateral displacement oscillation ^45-^45 and 
of the static stability for variations in the position of the towing 

point (1) in horizontal direction from the center of gravity toward 
the nose of the fuselage, and (2) in vertical direction from the center 
of gravity toward the wing nose. 
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Figure 10.- Period of oscillation and excitation time t* of 

45 45 

the lateral displacement oscillation and damping time of the static 
stability t- in the case of variation of the cable length. G = 5500 kg, 
6 
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Figure 12.- Period of oscillation and excitation time of the lateral 
displacement oscillation and damping time of the static stability- 
in the case of variation of the flying weight and for various flight 
speeds. k s = 6. 
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Figure 13.- Period of oscillation and excitation time of the lateral 
displacement oscillation and damping time of the static stability- 
in the case of variation of the climbing angle. G = 5500 kg, 
c a = 0.681, k s = 6. 
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Figure 14.- Roots of the lateral displacement oscillation U45 t UV45 
and of the static stability ^g in the case of variation of the dihedral. 
G = 5500 kg, c a = 0.436, k g = 6. 
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Figure 18.- Roots of the lateral displacement oscillation u^_g t iw^g 

and of the static stability in the case of forward deviation of 

the towing point at the lower side of the fuselage, a = horizontal 
component of the distance from towing point to the center of 
gravity. 






